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We examine the collective scattering of a superfluid droplet impinging on a mesoscopic Bose-Einstein con- 
densate (BEC) as a target. The BEC consists of an atomic gas with two internal electronic states, each of which 
is trapped by a finite-depth external potential. An off-resonant optical laser field provides a localized coupling 
between the BEC components in the trapping region. This mesoscopic scenario matches the microscopic setup 
for Feshbach scattering of two particles, when a bound state of one sub-manifold is embedded in the scattering 
continuum of the other sub-manifold. Within the mean-field picture, we obtain resonant scattering phase shifts 
from a linear response theory in agreement with an exact numerical solution of the real time scattering pro- 
cess and simple analytical approximations thereof. We find an energy-dependent transmission coefficient that is 
controllable via the optical field between and 100%. 
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The natural way to investigate quantum objects is scatter- 
ing. By selecting a convenient physical "stencil" with a large 
interaction cross section, one can probe the structure as well 
as the excitation properties of a target. It mostly happens that 
quantum objects are either microscopically small, like atoms 
or nuclei, or they are embedded inside a solid state system like 
electrons, electronic holes or Cooper pairs. Thus, they have an 
elusive character, which usually shuns direct observation. 

The generic response of a compound quantum object to 
bombardment with projectiles is either individual, i. e., by 
instantaneously ejecting another single particle of the com- 
pound, or it is collective, when after some transient period 
the target responds as a whole. Probably, the most drastic 
instance of this collective behavior is nuclear fission, when 
a heavy nucleus breaks up into large fragments - unleashing 
large amounts of kinetic energy. Modeling the dynamics of 
the atomic nucleus as a classical liquid drop OJEI gave a very 
intuitive interpretation of the observed phenomenon. 

Superconductivity in metals is another prominent collec- 
tive effect. Within the Ginzburg-Landau theory jj], one as- 
sociates a collective wave function with the Bose-condensed 
electronic Cooper-pairs, and a hydrodynamic description is 
again successful. The quantum mechanical nature of the fluid- 
like order parameter is usually discussed with the Josephson 
effect, but Andreev-Saint-James reflection ]4[ is an equally 
interesting phenomenon and much more in line with collec- 
tive scattering theory that will be presented in the following. 
This effect explains the unusual electrical transport properties 
through a normal metal-superconductor (N-S) junction. The 
occurrence of collective quantum mechanical resonances can 
be explained by a conversion of normal conductor electrons 
into hole-like excitations at the interface 0] ■ 

The discovery of superfluidity in bosonic 0] and, most 
recently, fermionic atomic gases MUM was an amazing 
achievement and is another manifestation of collective many- 
particle physics. Probably, the use of interatomic Feshbach 
scattering resonances has been the most fruitful novel concept 



of the past few years. Today, they are the universal tool to 

manipulate the binary interaction in an atomic gas 

in real time and they paved the way to fermionic superfluidity 

In the present Letter, we will demonstrate that the micro- 
scopic physics of binary Feshbach resonances can be also im- 
plemented at the mesoscopic level of an atomic BEC, giving 
rise to collective Feshbach resonances as shown in Fig.^ In 
particular, we will study the scattering properties of a weak 
coherent perturbation on a two-component BEC confined in 
quasi one-dimensional square-well potentials of finite-depth. 
Based on a Gross-Pitaevskii (GP) mean-field picture, we de- 
rive a two-component linear response Bogoliubov theory ifl^l 
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FIG. 1: Collective Feshbach resonance: complete transmission (A) 
or total refection (B) of a small coherent atomic wave-packet %pb (x, t) 
with incident momentum ko = 2.05 (A) or ko = 2.35 (B), when 
scattering-off a stationary, two-component BEC, trapped around 
— 1 < x < 1, with a maximal density 7ib(x) ~ 500, which is off- 
scale. The dimensionless density nt,(x, t) is depicted versus position 
x (in natural units of the trap lol i for three instants t: initially (ti, 
dotted), on impact (to, solid) and finally (tj, dashed). 
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for the scattering phases of the continuum perturbations. This 
is compared with the numerical simulation of the nonlinear 
wave-packet propagation in real time. Finally, we can ex- 
plain the appearance of collective Feshbach resonances qual- 
itatively from a simple Thomas-Fermi approximation of the 
Bogoliubov excitations. 

The dynamical response of a BEC has been the subject of 
intensive investigations, during the last decade. So far, laser 
light has been mostly the method of choice to impart momen- 
tum onto a BEC, to excite collective modes and to measure 
the dynamic structure factor fl7ll . However, light has also 
been used indirectly to prepare colliding matter-wave packets 
that have exhibited stimulated amplification as a result of their 
non-classical bosonic nature Hil . While the collision energy 
of these wave packets was rather low in this setup (few times 
the speed of sound), more recently, high energy, but ultra-cold 
thermal 17^1 and condensed atomic clouds 12011 of equal size 
have been used to measure the single-particle collision cross 
section. Theoretically, scattering of identical particles off a 
scalar BEC was already considered in the limits of low and 
hi gh e nergies I21H22I . Scattering of vortices was investigated 
in I23ll (already relating it to Andreev-Saint-James reflections 
01) and a generalization of Levinson's theorem was presented 
in El. 

Binary Feshbach scattering resonances cannot be described 
by a single channel potential scattering but they require at 
least the interaction of two different state manifolds. This 
was discussed first in nuclear scattering theory by H. Fesh- 
bach I25H26H27I1 and in the context of optical spectra of two- 
electron atoms by U. Fano l28l E^l . Thus, we will assume 
that the quasi one-dimensional condensate consists of atoms 
with two distinct internal electronic states. In the condensed 
phase, the system is described with order parameters ip a and 
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FIG. 2: Schematic set-up of a trapped two-component BEC with 
square-well potentials V a {x), Vb(x) (solid) and a coupling laser 
beam Sl(x) (long dashed) versus position x flUl . Scattering occurs 
in the two open (left/right) b-channels, while the a-channels are en- 
ergetically closed with a threshold energy A. The chemical potential 
is fi (heavy solid line), the dashed-dotted lines show the bound ex- 
citations and the numbered quasi-bound energy levels (dashed) are 
responsible for the collective Feshbach resonances. 



ipb that are coupled by an optical laser beam fl(x) with a large 
detuning A. Today, this can be achieved experimentally with 
very prolate traps, which freeze out the transverse degrees of 
motion, effectively. The corresponding two-component GP 
equation reads is also known from the interal Josephson effect 



idt + \d 2 x - 



n*(x) v™(x) 

V a (x)+g aa \lpa(x)\ 2 

V b (x) +g a b \4>a{x)\ 2 



gab \ipb(x)\ 2 
9bb \ipb(x)\ 2 



= 0, (1) 



where V a and Vb are external trapping potentials. For definite- 
ness, we pick the square-well potentials as they lead to simple, 
analytically solvable approximations. This setup is depicted 
in Fig. |2 The coupling constants g aa , g a b = gba, and gbb are 
proportional to the self- and cross-component s-wave scatter- 
ing lengths. For typical experimental values of 87 Rb see l3lll . 
but from the theoretical point of view the choice of parame- 
ters is uncritical 11511 . In order to describe the scattering of 
a superfluid droplet from the equilibrium BEC, we will de- 
termine the linear response modes of the two-component sys- 
tem 11611 . The stationary Bogoliubov ansatz with particle-like 
u = (u a , Ub) T and hole-like v = (v a , vt) T excitations is 

if>(x,t) = e-*"*^ ^) + e- l£t u(e,x) + e tet v* (e, x)}, 

where /.t is the chemical potential of the spinorial ground state 
xp^ = ( - 0a°' ) , tpb^ ) T and e is the excitation energy. This per- 
turbation to Eq. yields the four-dimensional Bogoliubov 
equations 
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The linear response energy matrix of Eq. (0 has the same 
structure as the well-known Bogoliubov equations in the case 
of a one-component condensate. It is symplectic and has a 
real-valued spectrum of pairwise positive and negative eigen- 
values I32I1 . Due to the finite-depth of the trapping potentials, 
the spectrum supports only a finite number of bound states and 
has a scattering continuum above a certain excitation energy. 
This is now the analogous situation as required for the two- 
particle Feshbach scattering resonances. If quasi-bound Bo- 
goliubov modes coincide in energy with continuum modes, 
we will obtain a resonance behavior. These bound, positive 
energy modes are depicted schematically in Fig.|2] For posi- 
tive energy solutions e > 0, these resonances appear in the do- 
main < e + /i = fc 2 /2 = E < A. In this regime, the mode 
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components u a (e, x), v a (e,x) and Vf,(e, x) are localized on 
the condensate and vanish exponentially for |x| — > oo. Only 
the particle-like component u&(e, x) can propagate outside 



lim Ub(e,x) 

; — >±oc 



cos(kx ± 8 e ), or sin(fca; ± S ). (3) 



Due to the reflection symmetry of the trapping potentials, 
the excitation can also be characterized by parity and we 
define even and odd phase shifts 5 e and S , respectively. 
These scattering phases have been evaluated numerically from 
Eqs. Jll2l3i and are displayed in Fig. [5] The analogy to the 
phenomenon of Feshbach resonances stands out very clearly. 
According to the Breit-Wigner parameterization of an isolated 
resonance 1 26], one finds approximately 



S les (E) = arctan 



r/2 



En — E 



2 
f 



dS tt 



dE 



(4) 



E=E r , 



In turn, one can determine the resonance energy Er and width 
r analytically from the poles of the S'-matrix in the vicinity of 
a Feshbach resonance 1331134113511. 




FIG. 3: Collective Feshbach resonances in the scattering phases 8 e 
(solid) and S (dashed) of a weak perturbation of a two-component 
BEC versus energy E. Note that the 7r -jumps of the collective Fes- 
hbach resonances occur at excitation energy Ei — fi + tj of the 
even/odd quasi-bound Bogoliubov modes of Fig. [2] The inset shows 
the Lorentzian behavior of the phase derivative S' a (E) close to the 
resonance as in Eq. J4j. For parameters see lol . 



From the phases shifts, one can obtain all scattering in- 
formation, like reflection R(E) or transmission amplitudes 
T(E), by considering a causal wave w[ (e, x), propagating 
to the right (k > 0), 



lim uf'(e,x) = e lkx + Re 



e l " cos(kx — 5 e ) + ie % " sm{kx — S a ), 



lim (e, x) = T e 



Ikx 



(5) 



e 15 " cos{kx + 5 e ) + ie 16 " sin(fcx + 5 Q ) 



Mr, 



The real-valuedness of the scattering phases of Fig.|3]implies 
a current conservation \R(E)\ 2 + \T(E)\ 2 = 1, and we can 
write the transmission coefficient in terms of the phases 



\T(E)\ 2 = co8 2 [6 e {E) - 6 {E)] 



(6) 



This transmission coefficient is shown in Fig.0] In the vicin- 
ity of the resonance energies E ps fj, + e^, it changes rapidly 
between and 100%. 
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FIG. 4: Transmission coefficient of a weak coherent perturbation in 
the b-component of the trapped BEC versus energy E from the lin- 
ear response Bogoliubov calculation (solid). The dashed line is the 
transmission coefficient obtained from propagating a Gaussian wave- 
packet in the nonlinear GP Eq. in real time. The dotted line rep- 
resents a crude Thomas-Fermi approximation. Feshbach resonances 
are seen clearly in all three curves superimposed on a background 
caused by potential scattering. 

In order to isolate the essential physical mechanism respon- 
sible for the collective resonance behavior, we approximate 
the GP Eq. Q in the Thomas-Fermi limit. Then ip^ (x) is 
constant within the square-wells and vanishes exactly else- 
where. In an additional approximation, we disregard the ma- 
trix M in the Bogoliubov self-energy Eq. (0. In this limit, the 
particle and hole excitations are decoupled and u satisfies a 
two-component Schrodinger equation 



(e + f-V B )u = 0. 



(7) 



Due to the TF approximation, this equation has again a square- 
well character and the solution can be found analytically 
i27i l33ll . Sparing the details of the calculation of the scat- 
tering phases, we simply present the results in the dotted line 
in Fig. |3 A good qualitative agreement needs to be acknowl- 
edged, while there are obviously shifts in the resonance ener- 
gies that are not accounted for in this simple approximation 
scheme. 

In a final step of the analysis of the collective Feshbach res- 
onance, we have also performed a numerical simulation of the 
nonlinear, time-dependent GP Eq. Q. We propagated an in- 
cident traveling Gaussian wave packet on top of the stationary 
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BEC solution depicted in Fig. [2 

M'J = 0) = + tf^e-^"*'. (8) 

It can be seen in Fig. ^that the b-component of the ground 
state solution is well localized in the trap center and has 
initially no overlap with the weak Gaussian perturbation if 
SN <C N and xq <C — 1. The initial momentum fco > 
of the wave was varied to cover the energy range in Fig. |4] 
In order to resolve the resonance structure, one needs a small 
momentum spread (a = 20) and we find that the wave-packet 
transmission spectrum matches the linear response approach 
very well. We show two instances of the propagation of an 
initial wave packet. In Fig. |T]\), the momentum fco « 2.05 
corresponds to an energy below the resonance energy, marked 
as Ei in Fig.|4] which leads to full transmission. In contrast, 
Fig. shows the total reflection of the wave-packet if the 
incident momentum fco ~ 2.35, corresponds exactly to the 
resonance energy Ei. 

In conclusion, we have identified collective Feshbach res- 
onances in a trapped two component BEC. These resonances 
do not require a binary Feshbach resonance to modify the bi- 
nary interaction, but quasi-bound Bogoliubov modes in the 
linear response spectrum. We have shown that this can be 
achieved easily with help of an optical laser beam. In con- 
trast to the microscopic binary Feshbach resonance, this quan- 
tum mechanical phenomenon exists on a mesoscopic, possibly 
macroscopic scale and can be used to control the transmission 
of matter-waves between and 100%. For conceptual sim- 
plicity, we have investigated a quasi one-dimensional geome- 
try with square-well potentials. None of this is important for 
the effect and it can be implemented in various experimen- 
tal configurations. However, studying this resonance in the 
mean-field picture can be only a first step and a proper inclu- 
sion of the thermal cloud is still lacking 136113711 . 
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